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Abstract 

We calculate the optical phonon correction to the binding energy of electrons 
to donors in cubic materials. Previous theories calculated the Rydberg energy 
reduced by the effective mass and the static dielectric function. They omitted 
an important energy term from the long-range polarization of the ionized 
donor, which vanishes for the neutral donor. They also omitted the donor- 
phonon interaction. Including these terms yields a new formula for the donor 
binding energy. 
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1 Introduction 



We present a calculation of the polaron correction to the binding energy Ed 
of electrons to donors. This topic was treated originally by Larsen[Tl IH [3l IH 
El E]. We adopt a similar approach, but include more terms in the theory. 

Prior theories are called the scaled hydrogen model (SHM). In SHM one 
takes the binding energy of hydrogen, the Rydberg energy, and scales it with 
the effective band mass m* of the electron, and the dielectric function e of 
the material. In SI units 

For weakly bound donors, the dielectric function is usually taken to be the 
zero frequency value b{0): which includes the polar contribution from op- 
tical phonons. When the donor binding energy is larger than the optical 
phonon energy, then some advocate using the high-frequency dielectric con- 
stant e{oo). Neither our theory, nor that of Larsen, were able to derive these 
expressions for the SHM. We show below that our theory gives the formula 
for the effective dielectric constant 
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The fractional factors are unexpected, and are the main result of our calcu- 
lation. The SHM is still in regular use today[71 El El [10]. 

In doing this calculation, we became aware that prior theories omitted 
an important energy term: the polarization energy of the ionized donor. An 
ionized donor is viewed as a point charge in a polarizable medium. There is 
a polarization energy associated with this long-range poential. We calculate 
below the exact expression for this energy. When the electron becomes bound 
to the donor, the donor becomes neutral, and the long-range polarization 
energy vanishes. It is replaced by short-range polarization on the scale of the 
Bohr radius of the donor. This difference in polarization energy contributes 
to the binding energy of the donor. It is one reason for the fractions in 
eqn.(2). 

We present calculations of donor binding energies for materials with the 
fee (face-centered-cubic) lattice. This includes most III-V and II- VI semi- 
conductors, and some oxides. For semiconductors with small donor binding 
energies, we get similar accuracy compared to the SHM. For the oxides, in 
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particular, our theory gives much larger binding energies than does the SHM, 
but neither agree with experiment. Our theory may also apply to the binding 
of holes to acceptors. However, the degeneracy of most valence bands makes 
this case more complicated, and we do not consider the acceptor case. 

Landau and Pekar[TT| [T2l [T3] developed the first theory of polarons for 
a free electron. They did a variational calculation on an electron eigenfunc- 
tion 0(r). They calculated the phonon displacement 5Q in response to the 
charge density of this eigenfunction, and then calculated the polaron energy 
oc {5QY. This phonon excess energy then became part of the variational en- 
ergy determining the electron eigenfunction. The Landau-Pekar theory later 
became recognized as the strong-coupling theory of polarons [H]. Here we 
adopt the same procedure to calculate the response of the phonons to the 
electron bound to the donor. Instead of the Gaussian eigenfunction used by 
Landau and Pekar, we use the exponential eigenfunction typical of a donor. 
Otherwise our calculation follows the same procedures. Obviously, our the- 
ory is also a strong coupling theory, which applies when £:(oo) and £:(0) are 
very different. 



2 Theory 

We employ an approach based on our earlier work|15t [TB] . The Hamiltonian 
includes optical phonons, donors at R; of charge g > 0, and an electron of 
charge e < 0, effective mass m*, which may be bound to a donor. The optical 
phonons have a displacement Qj, a Szigeti charge e*, which creates a dipole 
e*Qj-. 
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+ Vpd + Vep + Vpp (3) 

where: 

• The interaction between the phonons and the donors is 

_ e*q ^ Q, ■ (R, - RQ 
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The interaction between the electrons and the optical phonons is 
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where n(r) is the charge density of an electron bound to the donor. 
This reduces to the usual Frohlich interaction for electrons in a free 
energy band [17]. 



The dipole-dipole interaction between the phonons is 
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The lattice transform of this interaction is [15] 
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where Qq is the volume of the unit cell. We will usually take the limit 
of small wave vector in these expressions. 



3 Charged Donor 



First we evaluate the energy of an isolated donor, without an electron nearby. 
We transform the vibrations to collective coordinates 



The potential energy terms for LO phonons are 
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The second term in comes from the dipole-dipole interaction. The lattice 
transform of Vpd is 



qe 



47r£:(oo)VA^ k 



'tk-R, 



Ani k 



E^ R3 n 1,2 



(13) 
(14) 



where the latter limit is for long wavelength. In this case the expression for 
/3k is 



qe I 



^g.k.R, 



(15) 



The phonon displacements enter as {k ■ Qfc), where the notation empha- 
sizes that only longitudinal optical phonons are involved. The transverse 
(uto) and longitudinal (ulo) optical phonon frequencies are defined as 



07: 



TO 



LO 



Kt _ K 1 2 
Kl _K 2 ^ 2 



*\2 



(e*) 



MVLoe{oo) 



The static dielectric function is 
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We complete the square on the interaction in eqn. fllll) . 
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The last term is the effective interaction between the optical phonons and 
the donor, due to the phonon polarization 
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For 7^ R;/ the wave vector integral gives the effective interaction between 
the two donors 
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We combine this term with the direct interaction in eqn.([n]) to find 
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The latter identity follows directly from eqn. f[T8|) . The effective Coulomb in- 
teraction between ionized donors is screened by the static dielectric function. 

The case that / = I' gives the energy of a single donor from the phonon 
polarization. It has the approximate value of 
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where a is approximately a lattice constant. An exact expression for d is 
given below. This energy term comes from the polarization induced by a 
single charged donor, which reduces its energy. 



4 Electron Bound to Donor 

We repeat the above calculation but include an electron bound to a donor at 
R/. We do this calculation of the donor binding energy using a variational 
parameter a defined as 



^(^) = \ — 3exp[-ar/ao] (28) 



n{r) = |?/;(r)p (29) 
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where the Bohr radius is given by 



ao 
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and m* is the band effective mass of the conduction electron. Larsen[T] did 
a similar calculation but omitted the interaction Vpd which is the term that 
polarizes the lattice in the absence of an electron. The inclusion of this 
important term is crucial for including all essential interactions. Again we 
transform the vibrations to collective coordinates 



The potential energy terms in the limit of long wave-length are 
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The first term in 7^ comes from the donor-phonon interaction, and the second 
comes from the electron-phonon interaction. The donor-phonon interaction 
is evaluated in the long wavelength limit. We complete the square on the 
interaction, and also set q = —e 
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The last term is part of the effective interaction between the electron and 
the donor, due to the phonon polarization 

^ El7.r (36) 
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We evaluate this expression below in several limits. There are two types of 
terms. Ill ^ V then it is the phonon- induced interaction between two neutral 
donors. If / = /' it is the polaron correction to the binding energy of a single 
donor. Note that A(x = 0) = 1 so the integrand vanishes at /c = 0. There is 
no long-range Coulomb interaction between neutral donors. 
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4.1 Interactions between donors 



First consider the case that I ^ I' so the expression gives the phonon- 
induced interaction energy between two neutral donors. For the effective 
mass hmit, where the orbit radius oq covers many lattice constants, we can 
extend the limit of integration to infinity. Change to dimensionless variables 
X = kao/{2a),y = 2aRii'/ao, and we get 

2 2 

Vpp" = f , i{y) (38) 

Hy) ^ -I (39) 

The above integral can be evaluated by contour integration. It gives a term 

I{y) = —jZa^y^ (40) 

y m=0 

The interaction decays exponentially, which is typical of the interaction be- 
tween neutral charge distributions, when neglecting correlations. 



4.2 Neutral Donor #1 

We evaluate eqn. (l371) for the case that / = I'. Expand the factor 
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The first term ("1") we evaluate as eqn. (l271) . so the polarization energy of the 
donor is unchanged. The other two terms contribute to the donor binding 
energy. In the effective mass limit, the last two terms are evaluated assuming 
that the fc-integral goes to infinity. 
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The variational energy for the donor binding is 



£{a) = ERy[a^ - 2a{l - X)], X 
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(43) 
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Minimizing with respect to a gives 

f(a„) = -ER,(\-Xf (45) 

For weakly bound donors, the dielectric screening is not given by either Sao 
or e(0). Instead, it is 5/16 of £00 and 11/16 of £:(0). 

An interesting result is obtained if we neglect A^. In that case the total 
variational energy for the electron bound to the donor is 

^(«) = T^^-r^^a-^) (46) 



2m*al 47r£:(0)ao 
Varying a gives the minimum 



(47) 



«. ^ (^j ^ (48) 

The donor binding energy is screened by the zero frequency dielectric func- 
tion, which includes the contribution from the optical phonons. This is the 
usual form of the SHM, which is obtained only by neglecting A^. Since 
neglecting A^ is a poor approximation, the SHM is not accurate. 

4.3 Neutral Donor #2 

Here we evaluate the screening corrections to the neutral donor by a differ- 
ent method, which calculates all energy terms in real space. This has the 
advantage of not having to make a long-wavelenth approximation. The final 
expression should be more accurate. 
The interaction V^p is rewritten as 

Vep = . E ( E 1 E e"^-^-Qk • v^,/(i?,) (50) 
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When q = —e the 1/R term cancels Vpd- The combined terms give for 7^ 

PP* / \ ^ ( -2aR/ao ] 



The next step is to take X^k l7kP- There is a double sum over (j, j')- When 
/ = /' the summation over k forces Rj = Rj/. The final answer for / = /' is 

This answer is only slightly different than case #1. The summation over Rj 
converges to a finite value. There is no ionized donor energy as in eqn.f l27p . 
That is sensible, since the neutral donor should not attract long-range polar- 
ization. We rewrite the above expression using the nearest neighbor distance 
c? as a unit of length 

2 

V"pp = -ERy^^F{ad/ao) (53) 

F{ad/ao) = x^E p4 {I + 2aR, / + 2a' alf (54) 

An important result is the screening energy of the ionized donor, which is 
found by setting ao = 00 



This formula is exact, and is an improvement to eqn.f l27j) . For the fee lattice 
F{0) = 1.425780 using Ewald methods, which are described in the Appendix. 
By comparing the two expressions, we determine that 

/2 

a = a^-— = 0.992a (56) 

The energy of the ionized donor must be subtracted from the result of equa- 
tion ( 153|) . which gives the final donor binding energy 

£{a) = ERy{a' -2a-ri[F{ad/ao) - F{0)]} (57) 
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For a given crystal structure, F{x) is a function of x = ad/aQ. The self- 
consistent variational equation for the coupling constant at minimum 
energy is 

1 co^ 

"0 = 1 + 7;^r^'i^)^=0'od/ao (59) 

where F'{x) = dF/dx < 0. 

For many semiconductors the value of x is small, and it is adequate to 
take the limit of small x 

F{x)-F{0) = xF'(0) + O(x^) (60) 

Using Ewald methods we show for the fee lattice that F'{0) = —11/8, which 
makes eqn. fl57l) identical to eqn. fl43l) . The derivation in the appendix suggests 
that this result {F'{0) = —11/8) is valid for all lattices. For weakly bound 
donors our two derivations give the same result. Equations fl57|) and fl59l) are 
useful for materials with larger values of x. 

4.4 Atomic Limit for Donor 

One case is when the donor is tightly bound, so that d/a^ > 1. Even then the 
term F{x) is not small and cannot be neglected. For example, F{1) = 0.35. 

4.5 Ferroelectrics 

Many ferroelectrics have the feature that as one nears the transition temper- 
ature Tc then e{0) diverges while e{oo) remains constant. This is the result 
of uto going to zero, while ulq remains constant. In our first model, as 
Uto then Ulq ^i- The donor binding energy does not vanish, but is 
given by ^^^(5/16)2 ^ O.IE^^. 



5 fee Lattiee 

We present some calculations for the fee lattice, where d = a/\/2 and = 
a^/4 = d^/\/2. We discuss the evaluation of the function F(x), in eqn.( l5^ . 
and its derivative F\x). 

O fl p-'ixRj/d 

F{^) = '-pj:-^[^ + ^^R^/d + 2x'{R,/d)r (61) 
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Define the summation variable as Uj = Rj/d which starts out as one for the 
twelve nearest neighbors. 

F{x) = ^^^(l + 2x%+2x2y|f (62) 

4vrV2 yj 



One way to evaluate eqn. fl62|) is to sum all neighbors out to a distance y. 
However, the expression converges very slowly for small values of x. We 
performed the computer summations, out to a distance y^; and noticed that 
they seem to scale as —1/yx- We did a least squares fit to the expression 
F(0) = C — B/yrc, which gave C = 1.425781 which is very close to the Ewald 
result. Using Ewald methods in the Appendix we obtain the same value. 
A similar process determines F'{x) 

F\x) = --^;^^(l + 2xy, + 2x^2/1) (63) 

Using Ewald methods we found the exact expression -F'(O) = —11/8. In the 
appendix we derive the next terms in the power series in x 

F{x) = F{0)-—x + -x^ + O{x^) (64) 
8 3 

where B = 1.4595. With these values we constructed the graph of F{x) 
shown in figure 1. 

The eigenvalue equation (l59i) involves the function F'{x), which we have 
evaluated from equation (163|) . For fitting purposes, it is useful to have an 
analytic function that approximates these numberical values. The following 
expression is accurate to 0(0.3%) 

F'{x) ^ \e-^%l + 2x + 2x^) - e-^'^H + Bx^e'^^ 



7rV2 
11 



l+Cx + ^i^xf 



(65) 



where B is defined above, and the fitting parameter is = 5.09490. The 
first term on the right is the asymptotic limit, from the nearest neighbor 
contribution in eqn.f l63p . 

We have used our theory to calculate the donor binding energies in several 
crystals, as shown in table 1. Data is from Landolt-Bornstein[TS]. Only cubic 
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Figure 1: F{x) for the fee lattiee. The lower dashed curve shows the two 
first terms of the expansion in x, while the upper shows the results of the 
first four terms in x. 

crystals are included, and only those whose conduction band minimum is at 
the center of the Brillouin zone, so the effective mass is isotropic. For the ox- 
ides, the binding is sufficiently large that we solved eqn.( l59|) self-consistently. 
This gives a much larger binding energy compared to the SHM. Both our 
theory, and the SHM, are poor for the oxides, since electrons are too tightly 
bound for effective mass theory to be valid. We show these results only to 
demonstrate that for large values of x our theory is different than SHM. 

As stated in the introduction, we regard our theory as a strong-coupling 
result, which is valid when the two dielectric functions [e(oo) and ^(0)] are 
very different. This is not the case for most III-V semiconductors, and our 
theory does not do well for these materials. However, we expect our theory 
to apply quite well to the II- VI semiconductors, which are generally more 
polar than the HI- Vs. Table 1 shows our theory does relatively well for these 
materials. 

6 Discussion 

We have evaluated the polaron corrections to the binding energy of electrons 
to donors. We include the following interactions: electron-phonon, donor- 
phonon, and phonon-phonon. Two variational calculations gave identical 
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Table 1: Donor binding energies in meV. Only cubic crystals are included. 
Data from [TB] . 



results in the limit of a small binding energy. For the case of a large bind- 
ing energy, one method is exact, but gives an equation for the variational 
parameter that must be solved self-consistently. 

Most textbooks give that the donor binding energy as 

where m is the band effective mass of the electron, and £:(0) is the static 
dielectric function. We get a different expression 

e 16e(cx)) 16e(0) ^ ' 

which includes both dielectric constants [£:(oo), e(0)]. The effective Bohr 
radius is Oq = 47re jme^. We expect our theory applies in the strong- 
coupling limit, when £:(oo) and £:(0) are very different. 

There is also the question of whether the effective mass m* includes po- 
laron corrections, or is just the bare band mass. In our theory it is the bare 
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band mass. The above equation applies only to the case of Eo < Jk^lo- ^or 
materials with larger binding energies, one must solve a nonlinear equation 
to determine the parameters of the binding energy. We used our theory to 
calculate the donor binding energy of several materials with the fee crystal 
structure. 

One might consider a similar calculation for the electronic polaron effects: 
the terms that make e{oo) differ from Eq. Here one would start with a 
bare interaction, screened by the vacuum dielectric constant Eq, and consider 
how the electronic screening changes the donor binding energy. That is a 
different calculation than we have done here, since the electronic screening 
is in different places in the crystal. For the oxides, and other ionic crystals, 
the electronic polarization resides mostly with the anions[Tn]. However, for 
covalent materials, it resides in the bonds between ions. 

One of us (GDM) thanks the Erasmus Mundus program of the European 
Union for funding his trip to Sweden, where this project was initiated. He 
also thanks professors Per Hyldgaard and Elsebeth Schroder for hosting his 
visit to Chalmers University of Technology in Goteborg. 



Appendix: Ewald Summation 

The functions F{x), F'{x) are evaluated using an Ewald method [20]. We 
start with an evaluation of F{0) since it is easy. 



m = ^ 1:^ = 2^ EX (69) 

= Fl + Fr (70) 
An f^Jo ^ ^ 
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where rj = C/d. The term Fl we leave as is, since it converges rapidly in 
real space. The term Fr is changed to a summation over reciprocal lattice 
vectors by creating a periodic function of position r 

Vn{r) = r dtt'e-'^^'^^-^^' = Y.v{G)^''-' (73) 
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Table 2: Ewald summation for -F(O). 
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The first term in is the f (0), and the last term is from subtracting the 
Rj = term from Vr(0). Table 2 shows the separate contributions as a 
function of C . The final result is independent of C, which is a good way to 
check the computer code. 

Next we determine -F'(O) starting from eqnf l63p . Since the prefactor is 
proportional to x^, we must take the limit of x — )■ with some care. We use 
Ewald's method with an exponential 
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The lattice sum vanishes at x = : F[{0) = 0. Evaluate -Fr(x) by construct- 
ing a periodic function of r and then determining its Fourier coefficient. 

Vnir) = ^e-l''^-^l(^^'+*)[l + 2x|y,--r|+2x2(y,.-r)2] = ^e^G-'^;(G) 
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The result for -F'(O) comes from the term with G = 0. This integral is easy 
to evaluate 



v{0) = %Tiy/2 

r dtv(0) = 8nV2 
Jo 



1 Qx 24x2 



'j-S 



(83) 



-Qx' 



2 \16x2 (4x + r/)2 
1 1 ^ 



2x 



(4x)3 (4x + ?7)3^ 



(4a;)4 (4a; + r/)^ 



(84) 



The terms in the above equation, that are proportional to 0{x '^), give the 
limit as X — 
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One can also use the Ewald sum to find the first nonzero term in x which is 
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which gives B = 1.4595 ± 0.0001. 
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